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We reinvestigate the vaUdity of mapping the problem of two onsite interacting 
particles in a random potential onto an effective random matrix model. To 
this end we first study numerically how the non-interacting basis is coupled by 
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Q ' the interaction. Our results indicate that the typical coupling matrix element 



pg ' decreases significantly faster with increasing single-particle localization length 
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than is assumed in the random matrix model. We further show that even for 

^!^ ' models where the dependency of the coupling matrix element on the single- 
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t^^ ' particle localization length is correctly described by the corresponding random 

matrix model its predictions for the localization length can be qualitatively 

Vh ' incorrect. These results indicate that the mapping of an interacting random 

"A 

r^ ' system onto an effective random matrix model is potentially dangerous. We 
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J-^ ' also discuss how Imry's block-scaling picture for two interacting particles is 



influenced by the above arguments. 
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I. INTRODUCTION 

The interplay of disorder and many-body interactions in electronic systems has been 
studied intensively within the last two decades.lil For non-interacting electrons, the highly 
successful "scaling hypothesis of localization" was put forward in 1979 by Abrahams et 
alp but the role played by many-particle interactions is much less understood and still 
no entirely consistent picture exists.EJ The recent discovery of a metal-insulator transition 
in certain two-dimensional electron gases at zero magnetic fieldu has renewed the interest 
in this problem, since in the samples considered the electron interaction is estimated to 
be much larger than the Fermi energyB Thus the observed transition may be due to an 
interaction-driven enhancement of the conductivity. 

The simplest version of the interacting disordered particle problem is perhaps the case 
of just two interacting particles (TIP) in a random potential in one dimension (ID). For 
a Hubbard on-site interaction this problem has recently also attracted a lot of attention 
after Shepelyanskyoa argued that attractive as well as repulsive interactions between the 
two particles (bosons or fermions) lead to the formation of particle pairs whose localization 
length A, is much larger than the single-pattide (SP) localization length A J Based on 
a mapping of the TIP Hamiltonian onto an effective random matrix model (RMM) he 
predicted 

A2 ~ {u/vyxl (1) 

at two-particle energy E = 0, with V the nearest-neighbor transfer matrix element and U 
the Hubbard interaction strength. Shortly afterwards, ImryQ used a Thouless-type block- 
scaling picture (BSP) in support of this. The most surprising aspect of Eq. (0) is the fact 
that in the limit of weak disorder the ratio A2/A1 diverges. Thus, in the limit of weak 
disorder the particle pair can travel infinitely further than a SP. This should be contrasted 
with renormalization group studies of the ID Hubbard model at finite particle density which 
indicate that a repulsive onsite interaction leads to a strongly localized ground state.B 



Subsequent analytical investigations further explored the mapped TIP problem as an 
RMM problem.Bilj Direct numerical approaches to the TIP problem have been based on 
the time evolution of wave packets,a transfer matrix methods (TMM),t3 Green function 
approacheSjMlla or exact diagonalizationli^. In these investigations usually an enhancement 
of A2 compared to Ai has been found but the quantitative results differ both from the 
analytical prediction in Eq. (|l]), and from each other. Furthermore, a check of the functional 
dependence of A2 on Ai is numerically very expensive since it requires very large system sizes. 
Following the approach of Ref. |13|, two of us studied the TIP problem by a different TMME-3 
and found that (i) the enhancement A2/A1 decreases with increasing system size M, (ii) the 
behavior of A2 for f/ = is equal to Ai in the limit M -^ 00 only, and (iii) the enhancement 
A2/A1 also vanishes completely in this limit. Therefore it was concludedt3 that the TMM 
applied to the TIP problem in ID measures an enhancement of the localization length which 
is entirely due to the finiteness of the systems considered. 

In this paper we return the attention to the original mappinga of the TIP problem 
onto an effective RMM. We argue that the mapping as in Ref. ^ is potentially dangerous 
since (A) it overestimates the typical coupling matrix element and (B) it neglects phase 
correlations which we believe to be essential, because it is known that interference effects are 
responsible for Anderson localization to begin with. In order to establish that the mapping 
procedureP can lead to incorrect results we first numerically investigate the interaction- 
induced coupling matrix elements between the non-interacting basis states for various values 
of the SP localization length. We find that the typical coupling matrix element decreases 
significantly faster with increasing SP localization length than assumed in Ref. ^. This alone 
would lead to a significantly smaller increase, if any, of the TIP localization length than in 
Eq. (|1]). We further show that even if the RMM correctly described the dependency of the 
coupling matrix element on the SP localization length, its results for the TIP localization 
length cannot be trusted. To this end we present two simple physical examples, namely 
Anderson models with additional perturbing random potentials for which the RMM mapping 
yields the same enhancement of the localization length as for the TIP problem. However, 



for our examples this enhancement is obviously incorrect. We also show that analogous 
problems exist for the BSPE We argue that the failure of the RMM approach in our toy 
models is caused by neglecting the correlations between the coupling matrix elements. This 
has already been made responsiblet^'" for quantitative differences between Eq. (|1|) and 
numerical results for the TIP problem. We show, however, that neglecting the correlations 
not only changes the quantitative predictions of the theory but can lead to qualitatively 
incorrect results. 

The paper is organized as follows. In section || we briefly summarize the RMM approach 
to the TIP problem. In section |T| we present our numerical results for the TIP coupling 
matrix elements and their dependence on the SP localization length. The failure of the 
RMM approach to correctly predict the localization length of two toy models is discussed 
in section |I^ while section shows the failure of the BSP for these toy models. We discuss 



the relevance of our results for the original TIP problem and conclude in section VI. 



II. THE RANDOM MATRIX MODEL APPROACH 

Let us start by recalling the basic steps of the RMM approacha to TIP in a random 
potential. The relevant energy scales are chosen such that the SP band width AV is larger 
than the (uniform) spread of the disorder W which in turn is supposed to be larger than the 
interaction strength U. The basic idea is to represent the TIP Hamiltonian in the eigenbasis 
of the non-interacting problem and then to replace the full Hamiltonian by a suitably chosen 
random matrix. 

The (non-interacting) SP eigenstates are approximately described by 

(2) 



i{jn{x) ~ ^=^exp 
V Ai 



— + iOJx) 



Ai 
where Xn is the localization center of the nth eigenstate and On{x) is a phase which appears 

to be random but contains all the information about interferences necessary for Anderson 

localization. In the absence of interactions and neglecting symmetry considerations the 

two-particle eigenstates are just products of two SP eigenstates. 



tpnmix,y) ~ — exp 



+ i6n{x) +i9m{y) 



(3) 



Ai Ai 

where x and y are the coordinates of the first and second particle, respectively. Switching 
on the Hubbard interaction f/(x, y) = Ud^y between the two particles induces transitions 
between the eigenstates ipnm of the non-interacting problem. To estimate the transition rates 
it is first noted that the matrix element {ipnm\U\ipn'm') is exponentially small for \xn — ym\ > 
Ai or \xn' — ym'\ > Ai or |x„ — x„'| > Ai or \ym — ym'\ > Ai. Thus, the interaction 
couples each of the two-particle states (^ close to the diagonal in the 2D configuration 
space {\xn — ym\ < Ai) to 0{\l) other such states. The interaction matrix element is 
then the sum of Ai contributions each with magnitude f/A^^ and approximately random 
phases. Neglecting possible correlations among these contributions, Shepelyansky found the 
magnitude u of the matrix element 

Unmn'm' = (V'nm|f/|V^n'm') ~ UX^^''^, (4) 

independent of the interaction being attractive, repulsive or even random. Eq. (|^) is one of 
the essential ingredients of the RMM. In section |ITl| we will present numerical data in order 
to check its validity. We remark that the validity of Eq. (H) has recently been questioned in 



Ref. |T8| where the authors have computed a different estimate taking into account the nearly 
Bloch-like structure of the eigenstates for small W. 

ShepelyanskyQ'B now replaced the full TIP Hamiltonian by an effective RMM for those 
of the two-particle states that are coupled by the interaction. Thus the Hamiltonian matrix 
becomes a banded matrix whose elements are independent Gaussian random numbers with 
zero mean. The diagonal elements are drawn from a distribution of width V, because for 
small disorder W the nearest-neighbor transfer V determines the band width of the SP 
states. The distribution of the off-diagonal elements has width \U\\i within a band of 
width Ai. 

In order to obtain results for the localization properties of such an RMM one has to 
distinguish different regimes, depending on the strength of the interaction. If the interaction 



is strong enough to couple many non- interacting eigenstates, i.e., the inverse hfetime F of 
a non-interacting state is large compared to the level spacing of the coupled states, Fermi's 
golden rule can be applied. This regime was investigated in Ref. ^ and also gives the largest 
enhancement of A2 compared to Ai. We note that in this regime the level-spacing distribution 
of the non-interacting system cannot play a significant role since the interactions couple a 
large number of levels and lead to a decay into a quasi-continuum of final states. In the 
opposite limit, i.e., if the interaction couples only few non-interacting eigenstates, Fermi's 
golden rule cannot be applied. Instead, one finds Rabi oscillations between the few coupled 
states.tj We note that in this regime the level spacing distribution of the non-interacting 
states becomes important. In the following we will only consider the golden rule regime. 

The localization length of the effective RMM can be determined by several equivalent 
methods. Here we follow Shepelyansky:S Calculating the decay rate F of a non-interacting 
eigenstate by means of Fermi's golden rule gives F ~ U'^/XiV. Since the typical hopping 
distance is of the order of Ai the diffusion constant is Z^ ~ U'^Xi/V. Within a time r the 
particle pair visits A^ ~ U\^ y-1/27-1/2 g^^tes. Diffusion stops when the level spacing of the 
visited states is of the order of the frequency resolution 1/r. This determines the cut-off time 



T* and the corresponding pair-localization length is obtained as A2 ~ y/Dr* ~ (f//K)^Af 
in agreement with Eq. (|l]). Applicability of Fermi's golden rule requires F ^ V/Xf which 
is equivalent to U'^Xi/V'^ ^ 1. This is exactly the condition for an enhancement of A2 
compared to Ai. 

Let us recapitulate: The mapping of the TIP problem onto the RMM described above 
relies on two assumptions: (A) the non-interacting wavefunctions can be described by a 
decaying amplitude with finite localization length and a random phase which leads to the 
V jX-y behavior in Eq. (^) and (B) any correlations between the matrix elements in the 
Hamiltonian can be neglected. In the next two sections we will closer analyze the validity 
of these two assumptions. 



III. NUMERICAL RESULTS FOR THE MATRIX ELEMENTS 

In this section we present results for the interaction-induced couphng matrix elements 
in order to check whether they follow the Xi power law (^ as assumed in Ref. ^. Since 
Ai deviates from the simple power-law predictionO Ai ~ lOA/W'^ at E = already for 
■^1 ~ 4 {W > 5), we have first computed Ai by TMMpI in ID with 0.1% accuracy for 
all W > 0.3 (Ai ~ 1156). We next exactly diagonalize the SP Hamiltonian and obtain 
the eigenstates. We then compute the "center-of-mass" (CM) of these eigenstates as Xn = 
Xlx'^lV'nl^)!/ Xlx iV'nl^)!- For hard wall boundary conditions, we have checked that using 
this definition of the CM we can reproduce the disorder dependence of Ai from the decay of 
the SP wave function ipn via 1/Ai = — lim|a;_x„|^oo In |'?/'n(a^)|/|a; — Xn\ to within 10% up to 
Ai = 104 {W = 1) for 50 samples of length M = 1200. For periodic boundary conditions, 
we use a suitably generalized definition for the CM. We next calculate the matrix elements 
ii^nmlUlipn'm') for all statcs with appropriate CM, i.e., |x„ — ym\ < Ai, |x„/ — ym'\ < Ai, 
\xn — Xn'\ < Ai and \ym — Vm'l < Ai. Since the interaction strength U appears only as a 
multiplicative prefactor in the matrix elements, we choose U = 1 in all of what follows. We 
emphasize that the bottleneck in such a computation is not the system size M, but rather 
the exponentially growing number of overlapping matrix elements for increasing Ai. 

In Fig. |I| we show the unnormalized probability distributions Pd/o{u) of diagonal and 
off-diagonal coupling matrix elements. Pd/o{u) was computed at disorder W = 2 where the 
enhancement of A2 with respect to Ai is expected to be Iarge.ll2rt3 We have averaged over 
50 different disorder configurations for M = 200. For a more detailed inspection we plot the 
data on a doubly logarithmic scale in Figs. ^ and |^. As already discussed beforec3o'Ej we 
note that (i) the diagonal elements are non-negative, (ii) Po{u) is symmetric around u = 0. 
The deviation from symmetry for \u\ > 0.02, i.e., Pod^l) > Po{—\u\), is most likely due to 
the finite size of the samples. More importantly, (iii) Po{u) is strongly non-Gaussian. We 
remark that a fit to a Lorentzian distribution does also not describe the data, (iv) apart 
from a peak at m ~ 0, Pd{u) is approximately Gaussian, (v) Po{u) and Pd{u) have rather 



long tails, (vi) the total distribution of matrix elements P{u) is dominated by Po{u) (as in 
any matrix) and thus P{u) is strongly non-Gaussian with long tails. 

For such a non-Gaussian P{u) the average of the absolute matrix elements Ushs = {\u\), 
with (•) denoting the average over u according to P{u), is strongly influenced by rare events 
in the tails of the distribution. This is even more so when using the mean-square value 
^y(vF}. However, in the physical problem considered here these rare large couplings lead to 
oscillations of the system between the corresponding two TIP states but not to delocalization. 
The typical YdAne Utyp is thus better defined as the logarithmic average Wtyp = exp[(log(|M|))]. 

We have calculated both Uabs and Utyp for different values of W and 50 samples for 
M < 200 and 30 samples for M = 250. As shown in Fig. ^, the dependence of Wabs on Ai 
for Ai > 5 follows Wabs oc A^°. A fit for 20 < Ai < 111 yields a = —1.5 ± 0.1 as predicted in 
Ref. ^. However, the typical TIP matrix element Wtyp decreases much faster with increasing 
Ai. Fitting the Wtyp data to a power law for 20 < Ai < 111, we obtain a = 1.95 ± 0.10. 
Furthermore, for the largest Ai the numerical data deviate — albeit weakly — from the above 
power law showing a slight downward curvature in the double-logarithmic representation. 
In fact, if we only consider the data points from the large chains with M = 250, we already 
find a = 2.02 ± 0.10. This indicates that asymptotically the dependence is even stronger 
than A^^-^^ 

In Fig. 1^, we show Wabs and Utyp for the diagonal matrix elements only. For 20 < Ai < 111 
the data can be fitted by Md,abs oc ;\~o-9='=o-i and Md,typ oc A]^^'°''^°'^ Thus as expected Wd^abs 
and Md.typ behave similarly since Pd{u) may be approximated by a Gaussian distribution. 
Furthermore, a = 1 is in agreement with Eq. (21) of Ref. p!8| . 

Repeating the RMM calculation of section |I| with a dependence Utyp ~ U\i" instead of 
Eq. (^, we obtain 

A2 ~ (U/VfXf''^' . (5) 

If we use now use a = 1.95 ± 0.10 we find A2 ~ Xl'^^^"^. However, as discussed above the 
true asymptotic dependence of Mtyp on Ai is likely to be even stronger than A^^'^^ which 



in turn results in an even weaker enhancement of A2. We emphasize that the enhancement 
predicted by ShepelyanskyoS will vanish for a = 2 in the limit Ai -^ 00. A value of a > 2 
will in fact correspond to even stronger localization of the TIP. 

In order to further explore the validity of assumption (A) we compute Pd/o{u) for a 
site-dependent random onsite interaction U{x) G [—U,+U], averaging as before over 50 
samples. If assumption (A) is correct, the resulting distribution of the coupling matrix 
elements should qualitatively be similar to the one obtained for the original TIP problem. 
But as shown in Figs. |^ and |^, we find that the randomness of the interaction already leads 
to a significant decrease of the long-range nature of P{u). E.g., at \u\ = 0.02, there is a 
reduction in Pd/o(0.02) by a factor of approximately 10 for diagonal and approximately 5 
for off-diagonal matrix elements when compared to Pd/o(0.02) of the original TIP problem. 
We further compute P{u) for states with the same CM as previously, but otherwise chosen 
according to Eq. (|^) with uncorrelated random phases and exponentially decaying envelope. 
The disorder averaging is again over 50 samples. As shown in Fig. Q, Pd{u) now has a 
maximum at finite u. For these states Po{u) is well approximated by a Gaussian just 
as expected by Shepelyansky.Q'Q The double-logarithmic plot of Fig. |^ shows deviations 
from the symmetry Poiu) = Po{—u) for \u\ > 0.008, i.e., Po{\u\) > Po{—\u\). As for the 
TIP problem, we attribute this to the finite size of the samples considered. Again, we 
note that when compared to Po{u) for the TIP problem, the present distribution of matrix 
elements decreases much faster and at 1^1 = 0.02 is about one order of magnitude smaller. 
Furthermore, for |m| > 0.008 the distribution Po{u) is also smaller than that for the model 
with random interaction. Thus assumption (A) clearly oversimplifies the problem and the 
neglect of phase correlations leads to a wrong P^/o^u). In Fig. ||, we show -Uabs and Wtyp 
for the artificial states of Eq. (0). In complete agreement with our previous discussion, we 
find that for Ai > 20 both the average and the typical matrix element vary as m oc A^^'^ °'^ 
compatible with a = 3/2. 

In Fig. ^, we show also TIP data for chain lengths M = 100. We note that deviations due 
to the small system size lead to a smaller slope for Wabs and thus may give rise to an apparent 



enhancement of a. As can be seen in the figure, this decreasing of the slope happens for 
M = 100 aheady at Ai > 20 {W < 2.3). A power-law fit for 30 < Ai < 57 yields Uabs oc 
^-1.39 0.10^ rjj_^g finite-size deviations for Utyp are different. A power-law fit for 7 < Ai < 30 
gives Mtyp oc _x-i-'^'^±o-io whereas for 30 < Ai < 57 we find Utyp oc _x^2.06±o.io^ Thus first there 
is a decrease of a followed by a finite-size increase of a. For still larger Ai ^ M/2 the 
finite-size deviations of u^hs and Wtyp become very large even resulting in a positive slope. 
This finite-size effect may be at least partially responsible for the enhancement observed in 



Refs. |T3| and ^ for this value of M. 

In Fig. p, we show Uahs and Utyp for M = 100 with hard wall and periodic boundary 
conditions. Up to Ai ~ 10, the data for both boundary conditions agree quite well. For 
10 < Ai < 25, the slope for the data with periodic boundaries is slightly smaller than for the 
data with hard wall boundaries. Lastly, around Ai ~ M/2, the data for periodic boundaries 
shows a very fast decrease of u. Thus the data for periodic boundaries is influenced by the 
finite size of the sample already earlier than the data for hard wall boundaries. Nevertheless, 
except for these flnite size effects, our results for both boundary conditions are similar and 
we will restrict ourselves to the hard wall boundaries in the following. We remark that most 
numerical studies of the TIP problem also use this type of boundaries.li2rt3 

IV. FAILURE OF THE RMM APPROACH FOR TOY MODELS 

In this section we show that even an RMM which contains the correct dependence of 
the coupling matrix elements on the SP localization length may give qualitatively incorrect 
results. To this end we consider two toy models, viz. Anderson models of localization with 
additional perturbing random potentials. By a procedure analogous to that of section || we 
map these models onto RMMs and then show that these RMMs give erroneous enhancements 
of the localization length. 
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A. 2D Anderson model with perturbation on a line 

The first example is set up to lead to the same RMM as the TIP problem. It consists 
of the usual 2D Anderson model of localization perturbed by an additional weak random 
potential of strength U at the diagonal x = y in real space. Since this increases the width 
of the disorder distribution at the diagonal we expect the localization length to decrease. 
We now map onto an RMM as in Refs. ^J^. As above, the eigenstates of the unperturbed 
system are localized with a localization length Ai and approximately given by 



i^n{x,y) ~ — exp 



" ''-^ + ^9Jr) 



(6) 



Ai 

where r = (x, y)'^ is the coordinate vector of the particle and 6 is again a phase which is 
assumed to be random. The Hamiltonian of this model differs from the TIP Hamiltonian 
in two points: (i) the diagonal elements are independent random numbers instead of being 
partially correlated as in the TIP problem and (ii) the interaction potential U{x, x) G [—U, U] 
at each site of the diagonal is random instead of having a definite sign and modulus U as 
in the TIP problem. However, none of these points enters the mapping procedure outlined 
above. Thus, applying exactly the same arguments as for the TIP problem in section |I| 
we find that the perturbation couples each state close to the diagonal (|x„ — ?/„| < Ai) to 
0{Xl) other such states. The interaction matrix element is again a sum of 0{Xi) terms of 
magnitude U/Xf and random phases giving a typical value of UX^ . Consequently, our toy 
model is mapped onto exactly the same RMM as TIP in a random potential. 

As for the TIP case we now numerically check the relation between the coupling matrix 
element and the SP localization length Ai. We first note that the disorder dependence of 
Ai in the 2D Anderson model is no longer approximated by the simple power law cited in 
section |II[tJ In fact, Ai is usually much larger in the 2D case for the same value of W. Thus 
we compute estimates Ai(M) as a function of W for quasi-lD strips of finite strip width M 
with 1% accuracy by TMM. We remark that due to the self-averagingEl of l/Ai(M) this is 
equivalent to computing Ai(M) for many samples of M x M disordered squares. In Fig. 
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1^, we show data of Ai(M) as a function of W. We take Ai(50) to compute the couphng 
matrix elements. Since Ai(50) is always larger than for smaller system size, this choice 
only means that we sum over a few additional but very small terms when computing u. 
Next, we calculate both Wabs and Wtyp for different values of W and various M x M squares. 
Disorder averaging is over 20 samples and we study -Uabs and Wtyp as functions of Ai(M). We 
emphasize that instead of the well-known extrapolations of Ai(M) to infinite system size by 
means of finite-size scaling,c2l we take the finite-size approximants Ai(M) on purpose, since 
we compute A2 also for comparable finite sizes only. 

In Fig. 1^ we show the computed distributions Pd/o{u) for the present model. As for the 
TIP model the diagonal elements are non- negative and Pd{u) has a large peak at m = 0; 
Po{u) is again strongly non-Gaussian. The results for Wabs and Wtyp are presented in Fig. 
|. The dependence of Wabs on Ai(M) for 2 < Ai(M) < 12 follows Wabs oc Xi{M)~^-^^^-^ in 
agreement with our above prediction. Furthermore, here we also have Wtyp oc Ai(M)^^-^^°'^. 
As before, we note that the slopes of Uabs and Utyp become smaller for Ai(M) ^ M/2 due 
to the finite sample sizes. This finite-size effect is just the same as for TIP and thus further 
supports our use of the finite-size values Ai(M). We remark that if instead of Ai(M), we use 
Ai(50) for plotting the Wabs and Utyp data, that is irrespective of the system sizes for which 
they had been computed, we obtain Wabs oc _x-i-54±o.io ^^^ ^^^^ ^ j^-i.47±o.io^ Thus both 
choices of Ai show that Uabs and Wtyp vary as A]~^'^ within the accuracy of the calculation. 

Since our toy model is mapped onto the same RMM as the TIP problem the resulting 
localization length along the diagonal is also given by Eq. (p. We thus arrive at the sur- 
prising conclusion, that adding a weak random potential at the diagonal of a 2D Anderson 
model leads to an enormous enhancement of the localization length along this diagonal, 
in contradiction to the expectation expressed above, viz. that increasing disorder leads to 
stronger localization. 
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B. ID Anderson model with perturbation 

An even more striking contradiction can be obtained for a ID Anderson model of lo- 
calization. The eigenstates are again given by Eq. (j^) with Ai known from second order 
perturbation theoryEll and numerical calculationg3 to vary as Ai ~ V'^ /W"^ for small dis- 
order. We now add a weak random potential of strength U at all sites. Since the result 
is obviously a ID Anderson model with a slightly higher disorder strength the localization 
length will be reduced, Ai(f/) ~ V'^/{W'^ + f/^). Now we map onto an RMM according to 
Refs. ^^. The additional potential leads to transitions between the unperturbed eigenstates 
ipn- Each such state is now coupled to C(Ai) other states by coupling matrix elements 
{ipn\U\ipn') with magnitude u ~ f/A^ since we sum over Ai contributions with magnitude 
f//Ai and supposedly random phases. 

Again we numerically check the relation between Uabs and Utyp as functions of Ai. In Fig. 



10, we show results obtained for chains with various lengths and 50 disorder configurations 



for each W. Ai is computed by TMM as in section |T|. In Fig. ^ we show the distributions 
Pd/o{u)- We note that Po{u) is non-Gaussian as for the TIP model and the perturbed 2D 
Anderson model. Pd{u) is similar to the previous models, but the fluctuations are much 
larger. For 10 < Ai < 250, Uabs varies as ;\-o-48±o.io ^^ ^^ predicted above. Utyp varies as 
j^-o.59±o.io^ Both variations are compatible with a = 1/2. Again we need at least Ai > M/2 
in order to suppress the effects of the finite chain lengths. 

In analogy to section J| the application of Fermi's golden rule in this ID case leads to 
a diffusion constant D ~ U'^Xl/V. The number of states visited within a time r is now 
A^ ~ UXiV^^^'^r^^'^. Again, diffusion stops at a time r* when the level spacing of the states 
visited equals the frequency resolution. This gives r* ~ U'^Xf/V'^. The localization length A 



of the perturbed system thus reads A ~ \/ Dr* ~ t/^A^ as in Eq. (|l|), in clear contradiction 
to the correct result. 
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V. FAILURE OF THE BSP FOR TOY MODELS 

We now discuss the relation of our results to Imry's BSPQ for the TIP problem. In 
this approach one considers blocks of linear size Ai and calculates the dimensionless pair 
conductance on that scale, 

92 ~ ^, (7) 

where u represents the typical interaction-induced coupling matrix element between states 
in neighboring blocks and A ~ V/Xf is the level spacing within the block. If the typical 
coupling matrix element depends on Ai as n ~ UX^" the pair conductance obeys 

92 ~ {U/VrXl-'-. (8) 

Again, an estimate analogous to Shepelyansky's @ gives a = 3/2 which leads to a strong 
enhancement of the pair conductance g2 ~ Ai as compared to the SP conductance gi which 
is of order unity on scale Ai. In contrast, the numerical data of section |ITI| suggest that 
the pair conductance increases much less, viz. (72 ~ {U/VyXl'^^^'"^ for the fitted exponent 
a = 1.95 ± 0.10. Asymptotically for large Ai the pair conductance is likely to be enhanced 
even less than that. The behavior will be close to or even smaller than the marginal case 
(72 ~ 5'i- AH this is in complete agreement with our corresponding considerations for the 
RMM. 

For the 2D Anderson model considered in the last section, the BSP can be applied ana- 
loguously. Again, we consider blocks of linear size Ai and compute the typical perturbation- 
induced matrix elements between these blocks as in section |IVA] . We then find that accord- 
ing to the BSP the conductance of a 2D Anderson model with additional weak perturbing 
potential along the diagonal is given by Eq. (|^. Using a = 1.5 ± 0.1 as obtained in section 
p.V A| from the numerical data for Uabs and Utyp, we then have 92 ~ {U/V)'^Xi. Thus the BSP 
yields the same unphysical result as the RMM approach of section |IV A| . 

Let us also apply the BSP to the ID toy example. The level spacing in a ID block of 
size Ai is A ~ V/Xi, and the coupling matrix element between states in neighboring blocks 
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is t ~ UX^ . Thus, the conductance of the perturbed system on a scale Ai is obtained as 
dp ~ {U/V)'^Xi. For large Ai this again contradicts the correct result, viz. a decrease of the 
conductance compared to the unperturbed system. 

Thus, the BSP applied to the two toy models introduced in section |TV| gives the same 
qualitatively incorrect results for the localization properties as the RMM. This is not sur- 
prising since the only ingredients of the BSP are the intra-block level spacing A ~ V/Xl 
and the inter-block coupling matrix elements u which also enter the RMM and have been 
discussed in section |^ 

VI. CONCLUSIONS 

To summarize, we have reinvestigated the RMM approach to the problem of TIP in 
a random potential. We have shown that this kind of mapping an interacting disordered 
system onto an effective random matrix model is potentially dangerous since (A) it may 
overestimate the typical coupling matrix element and (B) it neglects correlations between 
the matrix elements. 

In the first part of the paper we investigated the dependence of the matrix elements 
entering the RMM on the SP localization length Ai. We found the dependence of the typical 
matrix element Utyp to be significantly stronger than for the averaged absolute value u^hs 
which is used in Refs. 3,|5|,[7|, |T0| , |l^ . If the RMM approach of section |l| is modified by using 
the numerically determined relation between utyp and Ai instead of Eq. (^) the resulting 
enhancement of A2 with respect to Ai becomes much weaker. We showed that the difference 
between Utyp and Wabs is due to the over-simplified assumption (A) that the wave functions 
behave according to Eq. (|^). Moreover, our data for Wtyp show systematic deviations from 
power-law behavior indicating that the true asymptotic dependence of Utyp on Ai is likely to 
be very close to or stronger than the marginal case Utyp ~ Xi^. If the asymptotic dependence 
is stronger than Utyp ~ A^^^ the Shepelyansky enhancement vanishes in the limit of large Ai 
even within the RMM approach. 
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In the second part of this paper we showed that there are physical situations where 
mapping onto an RMM as in Ref. ^ gives quahtatively incorrect results, e.g., an increase of 
the localization length in physical situations where it should rather decrease. This failure 
occurs even if the RMM contains the correct dependence of Wtyp on Ai. This shows in 
contrast to assumption (B) that in general the correlations between the matrix elements 
cannot be neglected since they contain information essential for the interference leading 
to Anderson localization. Note that the approach of Ref. |18|, while correcting assumption 
(A), still includes a mapping onto an RMM and thus is plagued by the same problems as 
assumption (B). Analogously, in Ref. ^the decay rate F is calculated numerically, avoiding 
assumption (A). However, the formula A2/A1 ~ V\\/V employed in Ref. |2^ is also based on 
an assumption similar to (B). 

Let us comment on the relevance of this work for the original problem of TIP in a random 
potential. None of our results constitute, of course, a proof that the enhancement of the 
TIP localization length A2 predicted in Ref. ^ does not exist. However, in our opinion, 
the toy counter examples to the RMM approach introduced in section |V| let the analytical 
arguments giving Eq. ([I|) appear much weaker. Taking the RMM approach seriously but 
using the numerical results for the typical coupling matrix element presented in section |T| 
we find that the dependence of the enhancement factor A2/A1 on Ai is significantly weaker 
than in Eq. (||). Nevertheless, an enhancement of the pair localization length for TIP as 
compared to Ai may still exist, although the underlying mechanism should then be different. 
Results supporting such an enhancement have been obtained by Green function methodalj 
together with finite-size scaling arguments.cJ The most recent data obtained in Ref. |15| finds 
an exponent a = 1.45 ± 0.2 for U = 1. 
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FIG. 1. Unnormalized distribution -Pd/o(^) of the diagonal (left panel) and off-diagonal (right 
panel) coupling matrix elements u with bin width A = 0.0003 for Ai = 26 iyV = 2) and M = 200. 
Circles and solid lines indicate TIP data, triangles and dashed lines indicate matrix elements 
computed using Eq. (0). The symbols mark the data for the 20 smallest \u\. 
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FIG. 2. Double-logarithmic plot of the unnormalized distribution Pd{u) of the diagonal cou- 
pling matrix elements u for Ai = 26 {W = 2) and M = 200 as in Fig. ffl. Solid, short-dashed, and 
long-dashed lines correspond to the TIP problem, the TIP problem with random interaction, and 
Eq. (g[), respectively. Circles, diamonds and triangles mark the data for the 10 smallest u in each 
case, only three of which are larger than 10^^ for the triangles. The vertical lines on the w-axis 
indicate Uabs (right) and utyp (left) as computed from the total distribution P{u). 
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FIG. 3. Double-logarithmic plot of the unnormalized distribution Po(^) of the off-diagonal 
coupling matrix elements u for Ai = 26 {W = 2) and M = 200 as in Fig. |I|. Solid, short-dashed, 
and long-dashed lines correspond to the TIP problem, the TIP problem with random interaction, 
and Eq. @, respectively. Circles, diamonds and triangles mark the data for the 10 smallest \u\ in 
each case. The vertical lines on the n-axis indicate Uabs (right) and utyp (left) as computed from 
the total distribution P(u). 
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FIG. 4. Dependence of Uabs {^, A) and utyp (o, t) on Ai for the TIP eigenstates (open 
symbols) and states chosen according to Eq. (^) (filled symbols) for M = 200. The small (bold) 
symbols indicate Uabs ^-nd utyp for AI = 100 {M = 250). The solid lines represent the power laws 
Uahs ~ K^'^ and Utyp ~ A^;^'^^. 
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FIG. 5. Dependence of Ud,abs (l^; ^) and tid,typ (°i ▼)> *-e., for diagonal matrix elements 
only, on Ai for the TIP eigenstates (open symbols) and states chosen according to Eq. (|3|) (filled 
symbols) for M = 200. The bold symbols indicate Ud.abs and Ud.typ for M = 250. The solid lines 



represent the power laws nd,abs ~ -^i ' and u^ 



l,typ ~ ^1 
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FIG. 6. Dependence of Uabs i^, +) and Utyp i°, x) on Ai for the TIP eigenstates for hard wall 
(□,o) and periodic (+, x) boundary conditions and size M = 100. 



25 



100 



J- 10 




10 



15 20 25 



w 



FIG. 7. Dependence of Ai(M) on disorder W for the 2D Anderson model at E = for 
M = 10, 25, 30, 35 and 50 indicated by increasing symbol size. We use the M = 50 data, emphasized 
by the solid line, as finite-size estimate of Ai. 
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FIG. 8. Unnormalized distribution for the diagonal (left panel) and off-diagonal (right panel) 
coupling matrix elements u with bin width A = 0.0015 for the perturbed 2D Anderson model with 
Ai = 3.1 (W = 12) and M = 25. The circles indicate the 20 smallest u (largest Po{u)) for diagonal 
(off-diagonal) data. 



27 



10" 



10" 



10" 




10 15 20 25 



^i(M) 



FIG. 9. Dependence of Uabs (squares) and Utyp (circles) on Ai(M) for the perturbed 2D 
Anderson model and M = 10, 25, 30 and 35 indicated by increasing symbol size. The solid lines 
represent the power laws Uabs ~ ^i ' and Utyp ~ Aj~ ' . 
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FIG. 10. Dependence of u^hs (squares) and utyp (circles) on Ai for the perturbed ID Anderson 
model and M = 200, 300, 500 and 800 indicated by increasing symbol size. The solid lines represent 
the power laws u^hs ~ -^r ^^"^ ""typ ~ A|f ' . 
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FIG. 11. Unnormalized distribution for the diagonal (left panel) and off-diagonal (right panel) 
coupling matrix elements u with bin width A = 0.0015 for the perturbed ID Anderson model with 
Ai = 26 (W = 2) and M = 200. The circles indicate the 20 smallest u (largest Po{u)) for the 
diagonal (off-diagonal) data. 
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